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We study the use of atom interferometers as detectors for gravitational waves in the mHz - Hz 
frequency band, which is complementary to planned optical interferometers, such as laser inter- 
ferometer gravitational wave observatories (LIGOs) and the Laser Interferometer Space Antenna 
(LISA). We describe an optimized atomic gravitational wave interferometric sensor (AGIS), whose 
sensitivity is proportional to the baseline length to power of 5/2, as opposed to the linear scaling of a 
more conservative design. Technical challenges are briefly discussed, as is a table-top demonstrator 
AGIS that is presently under construction at Berkeley. We study a range of potential sources of 
gravitational waves visible to AGIS, including galactic and extra-galactic binaries. Based on the 
predicted shot noise limited performance, AGIS should be capable of detecting type la supernovae 
precursors within 500 pc, up to 200 years beforehand. An optimized detector may be capable of 
detecting waves from RX J0806. 3-1-1527. 



PACS numbers: 

I. INTRODUCTION 

The production and propagation of gravitational waves 
is a central prediction of the theory of General Relativity. 
Direct observation of such waves has been attempted us- 
ing resonant bar detectors, which would be mechanically 
excited by passing gravitational waves, and using laser 
interferometers such as the laser interferometer gravita- 
tional wave observatories (LIGOs), in which the gravi- 
tational wave modulates the apparent distance between 
mirrors. These efforts have yet to detect gravitational 
waves, and improved interferometers such as LIGO-II are 
presently under construction, while the Laser Interfer- 
ometer Space Antenna LISA interferometer is being de- 
veloped. These experiments hope to sense gravitational 
waves by focusing the search upon waves with lower fre- 
quencies (LISA), where many sources of gravitational 
waves have already been optically identified, and/or by 
increasing their overall sensitivity (LIGO-II). 

In this article, we study the application of atom in- 
terferometry to build an Atomic Gravitational Wave In- 
terferometric Sensor (AGIS) [1], see Fig. 1. AGIS works 
on a similar principle to LIGO, replacing the measure- 
ment of the distance between macroscopic mirrors with 
measurement of the distance between freely falling atoms. 
These atoms can approximate an inertial frame with high 
accuracy. In contrast, the suspension of LIGO's macro- 
scopic mirrors requires elaborate seismic isolation sys- 
tems whose performance limits the sensitivity of LIGO 
at Hz-band or lower frequencies. 

AGIS benefits from the fact that atoms are all alike, 
and have relatively few degrees of freedom, eliminating 
many sources of noise that are problems for LIGO. There 
is no radiation pressure noise in AGIS, because each atom 
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FIG. 1: AGIS: two atom interferometers located on Earth or 
in space are separated by a distance L and addressed by the 
same laser radiation. 



interacts with a fixed number of photons. Thermal noise 
associated with the degrees of freedom of macroscopic 
mirrors and their suspension are likewise absent. Radia- 
tion pressure and suspension thermal noise are the dom- 
inant low- frequency noise sources in Advanced LIGO. 
These considerations suggest that AGIS might be more 
sensitive to gravitational waves at frequencies were the 
above noise sources are dominant, generally below 100 
Hz, complementing LIGO and LISA. However, AGIS 
must contend with a higher level of shot noise, since the 
flux of atoms in AGIS is much lower than that of photons 
in LIGO. This can nevertheless be mitigated by having 
each atom coherently interact with many photons. Grav- 
ity gradient noise remains important to both AGIS and 
LIGO. 

Thus, while AGIS seems interesting in the light of the 
above advantages, there are substantial challenges. The 
purpose of this study is to quantify some of these chal- 
lenges, in particular the ones that concern the atomic 
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physics aspects. We begin with a brief discussion of atom 
interferometers, and their use as key components of two 
AGIS detectors. We present both a basic configuration, 
relying on reasonable extrapolations from technologies 
presently available, and an optimized AGIS where pa- 
rameters are chosen to give maximal sensitivity for a 
given size of the apparatus. This sensitivity scales fa- 
vorably with the size of the apparatus. We then turn 
to a brief survey of the potentially detectable sources 
of gravitational waves in and near the AGIS frequency 
band. We give an overview of the state of the art of 
atom interferometry and on technologies that will need 
to be developed for AGIS; we also discuss some of the 
systematics that must be overcome. We hope this study 
will help to elucidate the promise and challenge of AGIS. 
Our study is mainly concerned with ground-based AGIS, 
though much of it will apply to space-based detectors [2] 
just as well. 

II. ATOM INTERFEROMETERS AND AGIS 

Light-pulse atom interferometry [3] has been used in 
measurements of local gravity [4] , the fine-structure con- 
stant [5-7], gravity gradients [8], Newton's gravitational 
constant [9, 10], a terrestrial test of general relativity 
that is competitive with astrophysics [11, 12], and a test 
of the gravitational redshift with part-per-billion accu- 
racy [13]. The field has recently seen revolutionary ad- 
vances in technology, improving the already impressive 
sensitivity of classical setups[14-17]. These technologies 
may be developed further into a tool for the detection of 
gravitational waves. 

Figure 2 shows two basic atom interferometer configu- 
rations, the Mach-Zehnder and Ramsey-Borde. In each 
case, an atomic matter wave is initially split by interac- 
tion with a light pulse, which transfers the momentum of 
n photons with a probability near 1/2. Subsequent pulses 
are used to redirect and then recombine the trajectories. 
When the matter waves from both paths interfere, the 
probability of observing the atom at a given output is 
given by (1 -f cos0)/2, where (j> is the phase difference 
between matter waves in both arms. This phase differ- 
ence 4> = 4>F + 4'i contains a contribution due to the free 
evolution of the wave function, (jip — ^Sci/h, given by 
the classical action Sc\ as evaluated along the trajecto- 
ries (A denotes the difference between both arms and h 
the reduced Planck constant). Another contribution (pi 
is because whenever a photon is absorbed (emitted), its 
phase is added to (subtracted from) the one of the matter 
wave. 

For a Mach-Zehnder interferometer, which will be the 
basic building block of AGIS, the leading order phase is 
given by 

(/)MZ = nkgT^, (1) 

where n is the number of photons transferred in the 
beam splitters, k is the laser wavenumber, g the lo- 
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FIG. 2: Left: Mach-Zehnder atom interferometer; Right: 
Ramsey-Borde Interferometer. Two conjugate interferome- 
ters are formed by interfering either the upper or lower out- 
puts of the second beam splitter. Trajectories which do not 
interfere are not shown. 

cal gravitational acceleration, and T the pulse separa- 
tion time (Fig. 2). For a Ramsey-Borde interferometer, 
4>KB = Sn^LUrT ± nkgT^, where the plus and minus sign 
refer to the upper and lower interferometer, respectively 
(Fig. 2) and ujr = hk'^/{2M) is the recoil frequency, 
where M is the mass of the atom. The first term arises 
because of the kinetic energy of the atoms due to the mo- 
mentum transfer of the photons. This term is absent in 
Mach-Zehnder interferometers, in which both trajectories 
receive momentum from the photons at some time. Re- 
cent technological advances in atom interferometry such 
as multiphoton Bragg diffraction, the use of simultaneous 
conjugate interferometers, and large momentum transfer 
using optical lattices may prove critical to the successful 
development of AGIS, and are reviewed below. 

A. AGIS 

A basic AGIS setup previously discussed by Dimopou- 
los et al. [1] is shown in Fig. 1. Two atom interfer- 
ometers, separated by a distance L, are addressed by 
a common laser system. A passing gravitational wave 
with strain amplitude h will modulate their distance L, 
and thereby the differential phase of the atom interfer- 
ometers. The differential phase modulation will have an 
amplitude of [1] 

A0 = 2nA:/iisin2(wT/2), (2) 

where u is the angular frequency of the gravitational 
wave. The shot-noise limit for the sensitivity of AGIS 
is thus given by 

" 2nkLsm^{ujT/2)^' 

where "q is the average atom flux through the inter- 
ferometer. The wavenumber k is determined by the 
atomic species used in the interferometer. For Cesium, 
k — 27r/(852 nm). Detection of gravitational waves will 
require the atoms momentum to be coherently split by 
thousands of Hk. Such large momentum splittings have 
yet to be experimentally demonstrated, but may be at- 
tainable by the use of accelerating optical lattices, as we 
discuss later in part IV C. 
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TABLE I: Parameters for AGIS that have been assumed in the plots in the theory section of this paper. 



Parameter 


Symbol 


Basic 


Optimized 


Wavcnumbcr 


k 


27r/852 nm 


27r/852 nm 


Momentum transfer/ (Rfe) 


n 


1,000 


31,000 


Pulse separation time 


T 


3s 


lis 


Tube length 


LTahe 


1,000 m 


3,000 m 


Separation 


L 




1,200 m 


Atom throughput 


V 


lO^Vs 


3 X 10"/s 


Peak sensitivity 


hrms 


7 X 10-2°/VHz 


1.3 X 10"^^/\/Hz 


Low freq. sensitivity 


LLF,opt 
' 'Tms 


3X 10-20(H5)2 1 


1.1 X 10-23(12)2 1 



AGIS will also require high atom throughput rj. 
Atomic fountains using Raman sideband cooling have 
demonstrated launches of 2.5 x 10^ state selected atoms 
at a three-dimensional temperature of 150 nK every two 
seconds [18]. As discussed in part IV A, we estimate 
that this flux can be scaled to provide throughputs of 
Tj = 10^^ atoms/s for AGIS, although shot noise limited 
detection of 10^^ atoms has not yet been demonstrated. 
We will use these parameters to define the characteristics 
of a "basic" AGIS detector, with a L = 1 km, as shown 
in Tab. I. This yields a root mean square sensitivity 
^rms = 7 X 10^^"/\/Hz at odd multiples of the frequency 
1/(2T) = 1/6 Hz, illustrated in Fig. 3. These parameters 
have been chosen because they appear feasible from ex- 
trapolation of the state of the art, at least in principle. 
Additional sources of systematic error are discussed in 
more detail in part IV E. 



-^Tube/2 [!]• The sensitivity is then 



^LF.opt 



4ff 



(5) 



or about 5 times better than with T = 3 s and the basic 
parameters in Tab. I. As discussed in the outlook, it may 
be possible to further extend the pulse separation time 
by trapping the atoms between beam splitting pulses. 

The size of the vacuum tube also limits the momentum 
transfer n, as the resulting spatial splitting nVrT, where 
Vr ^ 3.5mm/s is the recoil velocity of cesium atoms at a 
wavelength of 852 nm, must be accommodated. We con- 
servatively assume that the full spatial splitting must be 
added to i/xube- Globally optimizing the low frequency 
sensitivity with respect to T and n yields 



B. Optimized AGIS 

We can also consider strategics for optimizing this shot 
noise limit in free-fall configured AGIS, based on physical 
limitations. 

In the basic parameter set, we assumed an atom 
throughput of 10^^/s. As explained in part IV A, this 
may be increased to ~ 3 x 10^"^ /s using IkW of laser 
power for the two-dimensional magneto-optical trap (2D- 
MOT) . This higher flux can decrease shot noise by more 
than a factor of about 5. 

In a free-fall configuration, the maximum pulse sepa- 
ration time T is limited by the launch height of the atom 
interferometer. Optimizing T is crucial for reaching high 
sensitivity in the low- frequency limit w <C 27r/T, where 
the sensitivity is given by 



(4) 



The length Lxube of the vacuum tube for an AGIS with 
parameters L and T is about i-pube ~ i + gT'^/2, where 
g is the acceleration of free fall. Optimum sensitivity is 
reached when T — )■ -^Lxube/ff) when the launch height is 



^opt — 



«opt 
^LF,opt 




(6) 



This is 32 times better than the basic example of Tab. 
I for same tube length. Note the scaling of h^^f^^ with 

.5/2 
-^Tube- 

The length Lxube of the vacuum tube is limited mainly 
by the height of available mine shafts, or other facilities 
that can accommodate AGIS. 1 km has been assumed for 
the basic scenario in Tab. I. In principle, the deep under- 
ground science and engineering laboratory (DUSEL) at 
Lead, SD, is deep enough to accommodate a 3 km tube. 
Thus, Lxubc = 3 km has been assumed for the optimized 
scenario, giving another increase in sensitivity by a factor 
of 3^/2. 

All in all, the low-frequency shot noise limit of the 
"optimized" AGIS (Tab. I) is about 2,500 times as good 
as the basic one, a factor of 32 because of optimized n, T, 
a factor of 3^^^ due to Lxube) times a final factor of \/30 
due to the atom fiux. 
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III. SOURCES: BINARY INSPIRALS 

At least a third of all stars are in binary or higher 
multiple systems [19]. Binary systems are believed to 
be precursors to a variety of astrophysical phenomena of 
interest, ranging from type la supernovae [20], the forma- 
tion of neutron stars, and the evolution of supermassive 
black holes in galactic nuclei [21]. By detecting their 
emitted gravitational waves, AGIS might contribute to 
our understanding of such systems. The complete grav- 
itational wave spectrum emitted over the full lifespan of 
any given binary system is generally difficult to calculate, 
and particularly the final moments of the merger and 
ringdown phase when the system deviates significantly 
from the Newtonian approximation [22]. In contrast, the 
long Newtonian inspiral which leads up to the merger or 
collision of the two orbiting bodies is well understood, 
and can be adequately modeled with simple analytic ex- 
pressions applicable to a range of compact objects, in- 
cluding white dwarfs, neutron stars, and black holes. A 
binary system loses energy in the form of gravitational 
waves emitted with total luminosity [23] 



32 ^PM^ ^ , . 
-t^GW — -i=r-i_ E — J\£) 



(7) 



where /i — M1M2/M is the system's reduced mass, M = 
Ml + M2 the total mass, a is the distance separating 
the two stars, and G and c are the familiar gravitational 
constant and the speed of light, e is the eccentricity of 
the binary orbit, with /(e) given by [23] 



fie) = 



1 



73 

24^ 



37 
96 



(1-. 



2\-7/2 



(8) 



The angular frequency of the binary's orbit cj^ in the 
Newtonian limit, valid for all inspirals considered here, 
is — GM/a'^, allowing us to interchange a and ujb as 
needed. Binaries in circular orbits (e — 0) emit only at 
the second harmonic, but for non-circular orbits, (e 7^ 0), 
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FIG. 3: Root-mean-square strain noise per root Hertz for the 
two AGIS configurations described in the text and Tab. I. 



gravitational waves are also given off at higher harmonics 
WGW,n = nuJb, for n = 2, 3, 4, ... . The fraction F„(e) of 
the total luminosity Lq^^ emitted into the nth harmonic 
is given by [24] 



F„(e)=.g(n, £)//(£), 



(9) 



with 



gin,e) = — {[Jn-2ine) - 2eJ„_i(n£) + f J„(ne) 

+2eJ„+i(ne) - Jn+2{ne)f 
+{1 - [Jn-2{ne) - 2Jn{ne) + Jn+2{ne)f 



nine)]'}. (10) 



As the binary loses energy via emission of gravitational 
waves, the distance a between stars in a circular orbit 
decreases according to [23] 



a(i) = ao(l-^/ro)l/^ 



with 



To 



256 AiM2 



(11) 



(12) 



and ao — a{t = 0). This implies that the orbital fre- 
quency and gravitational wave luminosity of the system 
steadily increases in the Newtonian limit, unless other- 
wise perturbed (i.e. by collision or coalescence of the 
two bodies). Highly elliptical orbits emit more strongly, 
and inspiral more rapidly. Since most of the extra en- 
ergy is emitted at periastron, such orbits will gradually 
circularize over time [24, 25]. 

Gravitational waves can propagate with one of two or- 
thogonal polarizations, "-f" and "x". The power per 
solid angle emitted at a given polarization by a circular 
binary (e = 0) at an angle 6 to the binary rotation axis 
is 



^=^«-32 2;^(l + ^°^ ^) 
dPx , 5 2 2. 
ail 32 n 



(13) 
(14) 



The ability of a detector to see a source emitting one or 
both polarizations depends upon the orientation of the 
source relative to the distance vector to the detector, as 
well as the detector's geometry. For the purposes of esti- 
mating which objects AGIS might detect, we will assume 
that the detector is optimally oriented to detect the -I— 
polarized mode from any given source, and that the emit- 
ting binaries' rotation axes are rotated by 7r/3 from the 
distance vector linking them to the detector. The first 
assumption maximizes the signal in the detector, since 
> for all 9, while the second simplifies our anal- 
ysis, since the power emitted into the -I- mode is then 
equal to that which would be supplied by an isotropic 
and unpolarized source. 
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In the low energy density limit, where their energy 
density is insufficient to cause significant self-gravitation 
effects, gravitational waves propagate as solutions to the 
conventional wave equation. For a monochromatic wave 
with intensity $ and frequency w, the dimensionless am- 
plitude h of the wave is given by 

$ = ^(2Trufh^. (15) 

The gravitational wave is a disturbance of the metric of 
spacetime itself, acting tidally to alternately stretch and 
compress the distribution of matter and energy trans- 
versely to its direction of propagation [21]. For freely 
falling objects separated by a distance L, passage of an 
appropriately polarized gravitational wave with ampli- 
tude h causes the separation to vary sinusoidally with an 
amplitude of SL = hL/2. Interfcrometric detectors such 
as AGIS make precise measurements of such variations, 
and are thus sensitive to the amplitude of the gravita- 
tional wave, rather than its intensity. For the -I— polarized 
mode emitted at an angle 7r/3 from the quadrupole axis 
of the source, the amplitude at a distance r can be 
related to its total gravitational wave luminosity by 

$(r) = g(2..)^/.^ = i^, (16) 
so that the metric strain amplitude h is given by 

1 GLqW 

Note that since our detectors are sensitive to the strain 
amplitude h, rather than the energy density oc h'^, the de- 
tected signals drop as 1/r. A useful distance- independent 
definition of the "magnitude" of a particular harmonic of 
a gravitational wave source is thus most conveniently de- 
fined as rh [26]. 

To determine whether the magnitude of a given source 
of gravitational waves is sufficient for AGIS to detect, we 
convert our detector's /irms/v^Hz noise to a noise equiva- 
lent magnitude at the distance of the source. This noise 
equivalent magnitude at a distance r after an integration 
time T is given by 

r 

(r/l)nem(l') = -^/lrms('^)- (18) 

Note that detection may not be reliably accomplished 
until a significant amount of additional time has elapsed, 
depending upon what signal to noise ratio (SNR) thresh- 
old is required. For independent detection of gravita- 
tional waves against a Gaussian noise background, SNRs 
of 5 are typically considered necessary, while for signif- 
icantly non-Gaussian backgrounds, the SNR threshold 
may be higher. For joint-detection schemes, where some 
parameters of the emitted waves can be determined by 
e.g. optical observation, a somewhat lower SNR may be 
acceptable in some cases. Little is known about the noise 



backgrounds applicable to AGIS, so in what follows, we 
will simply plot the integrated noise equivalent magni- 
tude, and, assuming a Gaussian noise spectrum, require 
a SNR of 5 for detection. Objects within a distance r 
from the detector whose gravitational wave emissions lie 
above the plotted noise equivalent magnitude spectrum 
on the rh vs. p plot will be detectable with an SNR of 5 
after the observation period T. 

Detection of a gravitational wave requires continuous 
observations over at least one, and typically many, wave 
periods. The time required to accumulate a sufficient 
number of cycles to distingiiish the gravitational wave 
signal from the detector's background noise may in some 
cases exceed the lifetime of the source, precluding de- 
tection. To determine whether this is the case, we use 
Eq. (17) and Eq. (11), to relate the magnitude rh of a 
given £ = binary to the time r required to double its 
frequency. This relation is 

, , 75(8- 21/3) c ^ ^ 

(-^)<ioubie = 4096^7/2^^2 . (19) 

where v = 27ra; = 2cjb/ (27r) is the initial frequency of the 
emitted (second harmonic) waves. Note that the form 
of Eq. (11) is such that the time required for a binary 
to double its emission frequency is far greater than the 
lifetime remaining to it after its frequency has doubled. 
This c;haracteristic doubling time is independent of the 
constituents of the system damped by gravitational wave 
emission, and may therefore be broadly applied to the 
evolution of compact binaries composed of white dwarfs, 
neutron stars, or black holes. Using Eq. (19), we can 
easily determine whether a given source of gravitational 
waves will survive longer than the time required to ob- 
serve it. 



A. Compact Galactic Binaries 

Compact inspiraling binary systems located within our 
galaxy constitute a promising source of gravitational 
waves which may be detectable by AGIS. Such systems 
may involve white dwarfs (WD), neutron stars (NS), as 
well as black holes (BHs). We note that the discussion 
below assumes that AGIS will be shot-noise limited at 
frequencies near 0.01 Hz. This poses the challenge of 
overcoming gravity gradient noise, as will be discussed 
later. 



1. White Dwarf Binaries 

The number density of galactic white dwarf stars is 
estimated to be (4.8 x 10"-'') pc'^ [27, 28], of which 
at least a third may be expected to be involved in bi- 
nary or some higher multiplicity systems [19]. Although 
WD binaries are comparatively weak emitters of gravita- 
tional waves, their comparative ubiquity and proximity 
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to Earth makes them easily detectable by space-based 
detectors like LISA [21], and also potentially observable 
by ground-based AGIS detectors. Other ground-based 
detectors such as LIGO are primarily sensitive to signals 
above ~ 10 Hz. WD- WD and WD-NS binaries are in- 
sufhciently compact to exist at such frequencies, and are 
thus not detectable by LIGO. 

In contrast, an analysis [26, 29] of the gravitational 
wave strain magnitude rh of WD- WD and WD-NS bi- 
naries, shown in Fig. 4, reveals that AGIS detectors are 
sensitive to WD- WD binaries with period 200 s or below 
at a distance of 1 pc, and systems with period 27 s and 
below at a distance of 10 kpc. This is a range sufficient to 
detect most such fast binaries residing the nearer half of 
the galaxy. The blue-shaded region A in Fig. 4 shows the 
range of magnitudes and frequencies at which WD- WD 
and WD-NS binaries can emit, taking the Chandrasckhar 
mass to be Mch = 1.44Mq, while the green-shaded re- 
gion B indicates the magnitude of NS-NS binaries, where 
the neutron star masses are assumed to be between M^h 
and 2.4M0. 

A recent study [20] suggests that type la supernovae 
arising from WD- WD binaries with total mass greater 
than 1.4M0, where Mq is the mass of the sun, occur in 
our galaxy at a rate of 1 x 10~^ yr~^. Figure 4 shows 
that basic AGIS (see Tab. I and below) should be able 
to detect gravitational radiation from la supernova pre- 
cursors with SNR greater than 5 at ranges exceeding 500 
pc after five years of integration, and is sensitive to them 
as early as 200 years prior to the pending collision. 

AGIS may also be sensitive to gravitational waves at 
higher harmonics that are emitted from systems with 
even longer periods, provided they have sufficiently ec- 
centric orbits. Althoiigh most galactic WD binaries arc 
expected to be in nearly circular orbits, recent simula- 
tions [33] suggest that binaries residing in nearby glob- 
ular chistcrs may have higher eccentricities, due to in- 
teractions with nearby stars. There are at least 79 such 
globular clusters within 10 kpc of the Sun [33, 34]. Fur- 
ther study will be required to determine whether such 
globular clusters are likely to contain binaries with sig- 
nificant eccentricities in the AGIS detection band. 

Although the reference AGIS detector discussed above 
is not sufficiently sensitive to detect any presently known 
compact binaries, an optimized detector with the opti- 
mized parameters listed on Table I would be. Figure 5 
depicts five times the noise equivalent strain magnitude 
of the shot noise in such a detector at 1 kpc, relative to 
the magnitude of the "brightest" known binaries. This 
upgraded AGIS detector could detect gravitational waves 
from RX J0806.3-I-1527 with SNR greater than 5 in five 
years [31]. 



2. Neutron Star and Low to Intermediate Mass Black Hole 

Binaries 

NS-NS, NS-BH, and BH-BH binaries are attractive 
sources of gravitational waves that AGIS might study. 
These systems are ultra-compact, and can in many cases 
exhibit a simple slow inspiral over AGIS's entire detec- 
tion bandwidth. They are also massive, making them 
significantly easier to detect at a given distance as com- 
pared to systems involving WD stars. Unfortunately, 
they are also much less common. There are at present 
at least five known galactic NS-NS binaries that are ex- 
pected to merge within several hundred Myr: PSR J1756- 
2251, PSR J1906-F0746, PSR B1534-I-12, PSR B1913-16, 
and PSR J0737-3039 [35-38]. Of these five, only PSR 
J1906-F0746 and PSR J0737-3039 have sufficiently short 
periods to come within three decades of AGIS's peak sen- 
sitivity bandwidth. As illustrated in Fig. 6, the reference 
AGIS configuration considered here may be capable of 
detecting all NS-NS binaries with SNR 5 within 1 kpc 
that have less than 40 years before coalescence after one 
year's observations. The expected coalescence rate for 
our galaxy is not well known; one estimate is at about 
10-5 yr-i [39]. 

The plot in Fig. 6 of 5 times the noise equivalent mag- 
nitude at 16 kpc after a year's integration of the basic 
AGIS considered here shows that one year's observations 
may be sufficient to detect most galactic NS-NS bina- 
ries with total mass greater than AMq with SNR > 5 
at least a year before merger. Since 16 kpc/Vuionth ~ 
50 kpc/y^, the same curve indicates that binaries com- 
posed of a IOM0 black hole and a neutron star, or more 
massive partner, may be detectable more than two years 
from merger anywhere in the galaxy or Large Magellanic 
Cloud after one year of observations. 

Recently, an intermediate-mass black hole (IMBH) be- 
tween 500 and 10^ solar masses has been observed in the 
galaxy ESO 243-49 [40], about 92 Mpc from Earth. In- 
spirals of IOMq BH around this IMBH around this BH 
should be detectable by the optimized AGIS a month or 
more so before merger, as indicated in region H of Fig. 9. 
If the mass of the IMBH exceeds 2OOOM0, NS and WD 
inspirals may also be detectable. It must be noted, how- 
ever, that this particular BH is believed to be the stripped 
nucleus of a small dwarf galaxy, and may no longer have 
a significant number of compact partners. Its existence 
nevertheless suggests that there may be more IMBHs in 
the vicinity, making them promising candidate sources. 



B. Inspirals at the Galactic Core 

At the center of the galactic core, approximately 8 kpc 
distant from the Sun, there is believed to be a supermas- 
sive black hole with total mass M = 4 x 10*^ Mq [41]. A 
basic shot noise limited AGIS detector would be capable 

of detecting inspiral waves from compact objects in close 
orbit about the core, provided that such objects exist and 
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FIG. 4: White dwarf and neutron star binaries. Gravitational strain magnitudes of potential WD- WD binaries (blue-shaded 
region A) [26] and NS-NS binaries (green-shaded region B). Downward sloping solid lines indicate the AGIS noise equivalent 
strain magnitude spectrum at the indicated distances per root observation time for the basic detector parameters in Tab. 1. 
Objects with gravitational strain magnitudes larger than the equivalent noise magnitude hrms at a particular distance will lead 
to a signal to noise ratio (SNR) of 1 after the indicated period of integration. Sources which produce a SNR of 5 lie above 
the 5/irms curves and are potentially detectable by AGIS. Black circles on the boundary of the WD- WD binary region indicate 
points at which the components of WD- WD systems with total mass M come into contact with one another. Downward-sloping 
dashed lines represent the rh threshold above which inspiraling sources must double their emission frequency in less than 200 
years (lower line) or one year (upper line). Points in the lightly shaded region represent the gravitational strain magnitudes of 
the binary systems which make up the LISA verification sources, as given in [30] and [31], and also the more recently studied HM 
Cancri system [32]. These systems which include AM CVn systems (circles), detached WD binaries (squares), ultra-compact 
x-ray binaries (diamonds), and cataclysmic variables (triangles). Sources marked with circles are believed to be within 500 pc 
of the Sun. 



are close to merger. As indicated in Fig. 7, AGIS is sen- 
sitive with SNR = 5 to waves emitted from M = 2QMq 
black holes more than 3 years before coalescence with the 
galactic core, and to the waves emitted by a M = ZMq 
neutron star in its final year with slightly reduced SNR. 
Circular inspirals involving the galactic core will not emit 
at frequencies exceeding z^max, given by 



SGA/total 



2 mHz, 



(20) 



for A/total — 4 X 10^ Mq. Inspiral waves at higher fre- 
quencies are precluded as they would require the orbiting 
bodies to be within 2 gravitational radii of one another, 
where the Newtonian approximation is expected to be 
invalid. Inspirals from systems with lower total mass can 



form closer orbits, and thus emit at higher frequencies. 
Such low mass black holes are also expected to abound 
in the vicinity of the galactic core, with populations of 
the central parsec estimated to be of order lO"* [42-44]. 
Any such binaries which will merge in ~ 5 years may 
be detected at 8 kpc. Estimates for the capture rate for 
a given black hole span a disquietingly large range from 
5 X 10~^ yr"^ [45] and 10^^ yr^"'^ (and, for anomalously 
heavy neutron stars even 10~*yr~^) [46] due, in part, to 
the lack of realistic agreed-on models for the structure of 
galactic nuclei. This puts the detection probability for 
AGIS anywhere between order-unity and 10~^ per year. 
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FIG. 5: Enlargement from Fig. 4. 5 times the noise equiva- 
lent strain magnitude for an optimized L = 3 km AGIS de- 
tector at 1 kpc. The enlarged point indicates a known source, 
RX J0806. 3+1527, which would produce a signal with SNR 
greater than 5 after one year of integration. 

C. Extragalactic Extreme Mass Ratio Inspirals 

Intermediate mass ratio inspiral (IMRI) and extreme 
mass ratio inspiral (EMRI) waves from other galaxies in 
the Local Group could lead to signal to noise ratios of 
larger than 1 albeit not above the minimum of 5 that we 
assume is required for detection, as indicated in Fig. 8. 
Such sources might, however, be detectable at larger 
SNRs with the optimized AGIS. Elements of the Local 
Group lie at distances as low as 22 kpc for the nearer 
satellites of the Milky Way, out to distances in excess 
of 1 Mpc. There are two spiral galaxies other than the 
Milky Way in the Local Group, Triangulum (M33) and 
Andromeda (M31), both of which lie ^ 900 kpc away [47]. 
Using a detector with the optimized parameters listed in 
Table I would extend the range at which such systems 
could be detected to Gpc scales, potentially allowing us 
to study the gravitational wave spectrum of the Local 
Supercluster, as indicated in Fig. 9. It would also be 
capable of detecting signals from EMRIs involving the 
SMBH at the core of Andromeda, believed to have mass 
Af = 2 X WMq [48]. 

IV. TECHNOLOGY AND SYSTEMATICS 

We now turn to an analysis of the technologies required 
for the development of AGIS, and a discussion of various 
sources of systematic error. Detection of gravitational 
waves will require hundreds or even thousands of hk mo- 
mentum splitting. To date, large-momentum transfer 
has been accomplished via multiphoton Bragg diffrac- 
tion, discussed in IV B, to achieve smaller momentum 
splittings. Although present interferometers are very far 
from this level of performance, techniques involving ac- 
celerating optical lattices, such as the Bragg-Bloch-Bragg 
beam splitter, described in part IV C and in more de- 



tail in [17], may achieve the required splittings. AGIS 
will also require high atom throughput, the prospects for 
which are covered in part IV A. 

Other technical requirements include the minimization 
of wavefront distortions in the laser beams. This will re- 
quire a mode filtering cavity and high-quality optical el- 
ements located inside the vacuum chamber. Moreover, 
smooth frequency ramps for driving Bloch-oscillations 
are necessary, which can be generated by state of the 
art direct digital synthesizers. A fraction of the atoms 
may miss one or more of the momentum transfers and 
cause spurious interference, which must be eliminated in 
order to achieve shot noise limited signals. 



A. Atom Flux 

AGIS will also require high atom throughput. Atomic 
fountains using Raman sideband cooling have demon- 
strated launches of 2.5 x 10® state selected atoms at a 
three-dimensional temperature of 150 nK every two sec- 
onds by loading from a vapor cell MOT of about 7 x 10® 
atoms in a roughly 3mm-diameter cloud [18]. This flux 
can be increased by using a two-dimensional MOT, which 
typically achieve a flux of about 6 x 10^° atoms per sec- 
ond with a total laser power of about 0.6 W [49]. Scaling 
linearly, we can expect a flux of 10^'^ — lO^'^ atoms with 
a laser power of (10-1000) W, of which about 1/3 can be 
launched and cooled if the efficiency of Ref. [18] can be 
reproduced. Of course, such scaling may not be straight- 
forward: the laser power that can be achieved with a 
single commercial tapered amplifiers reaches about 2W, 
although 5 W tapered amplifiers are under development 
[50]. The 2D-M0T beams can be generated by many 
such chips (as single mode beams are not required for 2D 
MOT operation) , and the atom fiux of several 2D MOTs 
can be combined. 

Increasing the atomic density beyond 3 x lO^/cm"^, 
demonstrated by Treutlein et al, is undesirable because 
of mean field shifts [51], so a sample of 7 x 10^^ will have a 
10 cm diameter. Such a diameter is compatible with the 
thick beams required for AGIS in order to achieve a large 
Rayleigh range. Thus, while technical issues will need to 
be addressed, we are confident that an atom flux of more 
than f] = 10^^ atoms/s can be achieved for AGIS. 

Finally, atom detection methods with sufficient signal 
to noise ratios to detect 10"'^^ atoms at the shot noise limit 
need to be developed. This may involve high solid an- 
gle light collection optics for fluorescence detection, the 
use of charge coupled devices (CCD) cameras and image 
processing techniques to suppress stray light. Most im- 
portantly, laser frequency and intensity stabilization is 
required. Chopping or modulating the beam at frequen- 
cies above the 1/f noise floor may be used to overcome 
technical noise. Further sources of systematic error are 
discussed in more detail in part IV E. 
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FIG. 6: Neutron star binary sources. Detectability of NS-NS binaries (with masses between IAMq and 2.6Mq) and those 
involving black holes {M < lOM©). The upper shaded region C encompasses BH-BH and BH-NS binaries composed of a 
IOM0 BH with a partner whose mass lies between 1.44Mq and lOOM©. The three downward sloping dashed lines indicate the 
thresholds above which a source's orbital frequency doubles in less than 100 years (lower), two years (middle) and one year 
(upper). As PSR J0737-3039 lies approximately 600 pc away from Earth, GWs it emits will, in about 84 million years, exceed 
the AGIS's background noise after one year of observation, some 500 years before entering LIGO's detection band. 



B. Multiphoton Bragg difFraction 



Large momentum transfer by accelerated 
optical lattices 



Whereas classical atom interferometers use Raman 
transitions to transfer the momentum of two photons 
to the matter waves, multiphoton Bragg difFraction can 
be used to form large-momentum transfer beam split- 
ters. We have interfered cesium matter waves that were 
split by a momentum difference Ap of up to 24/ifc, the 
highest so far [14] (Other methods, also based on Bragg 
diffraction, have meanwhile achieved similar momentum 
splitting [15]). The sensitivity of interferometers rises 
proportional to Ap in measurements of inertial effects, 
e.g., local gravity, the gravity gradient, or gravitational 
waves. The sensitivity even rises proportional to (Ap)^ in 
other applications, e.g., measurements of the fine struc- 
ture constant a and the recoil frequency. 



With Bragg diffraction, increasing Ap beyond approxi- 
mately 24:hk is difficult, because the required laser power 
rises sharply with momentum transfer. We overcome this 
limitation by coherent acceleration of matter waves in op- 
tical lattices (Bloch oscillations, Ref. [52]) to add further 
momentum. With this Bloch-Bragg-Bloch (BBB) beam 
splitter, we have achieved Ap = 88/ifc; in simultaneous 
conjugate interferometers (SCIs) with Ap — 2Ahk, the 
BBB splitter allows us to see interferences with ^ 30% 
of the theoretical contrast, compared to ~ 2% with Bragg 
diffraction [17]. It is worth mentioning that the SCIs with 
BBB splitters use, all in all, 6 Bragg diffractions and 24 
optical lattices, see Fig. 10. Since the BBB splitter does 
not require higher laser power for increasing Ap, we ex- 
pect that technical improvements, discussed below, will 
allow us to reach a splitting of hundreds or even thou- 
sands of photon momenta. Gravitational wave detection 
aside, the BBB interferometer may also enable measure- 
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FIG. 7: Galactic core sources. The red-shaded region D at the center left represents the gravitational strain magnitudes of 
extreme mass capture inspirals which may occur in the galactic core. This involves a black hole with AIbu = 4 x W^Mq 
and in-falling compact objects with masses ranging from 0.5Mq (which form the region's lower bound) to lOOM© (forming 
the region's upper bound). The line which parallels the shaded region's upper boundary indicates the signal generated by a 
IO^Mq black hole in orbit about the larger black hole, while the line inside the shaded region indicates the signal generated 
by a similarly orbiting 3Mq neutron star. As before, the two dashed blue lines indicate the point at which inspiraling bodies 
double their frequency: in less than 5 years (lower dashed line) or in less than one year (upper dashed line). 



ments of the Lense-Thirring effect [53] , tests of the equiv- 
alence principle at sensitivities of up to Sg/g ~ 10"^'' 
[54] , atom neutrality [55] , or measurements of fundamen- 
tal constants with sensitivity to supersymmetry [56]. 

Other necessary advances include high-power, ultra- 
low noise lasers [57], increased interferometer areas (and 
hence sensitivity) [54], advanced algorithms for data 
analysis [58, 59], and new atom-optics tools. 



D. Simultaneous Conjugate Interferometers 

The sensitivity of atom interferometers is often limited 
by the effects of vibrations. Using SCIs has allowed us to 
cancel the influence of vibrations [16]. This allowed us to 
increase the time T between light pulses to 50 ms from 
1 ms for atom interferometers with Ap — 20/ifc, which 
corresponds to an increase in sensitivity, by a factor of 50 
for recoil frequency measurements and 2,500 for gravity 
measurements . 



E. Systematic Influences 

So far, only atom shot noise has been considered in this 
paper. Many other limiting influences have been studied 
[1], and further studies will be required. This, however, 
is beyond the scope of the present paper. 

The following paragraphs mention some of the require- 
ments for reaching the shot noise limit in the basic sce- 
nario of Tab. I. This shot noise limit may be expressed 
as a 1 firad phase uncertainty per root Hz in the inter- 
ference fringes of the atomic matter waves, or, at a cycle 
time of 2T — 6s, about 0.4 /irad in one cycle. For the op- 
timized scenario, these figures are 0.2 /irad per root Hz 
and 0.04 /irad per cycle, respectively. 

1. Magnetic fields 

Atoms in mp = quantum states only exhibit a 
quadratic Zeeman effect, 427 /iHz/(mG)^ for Cs. If we 
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FIG. 8: Signals from extragalactic sources. The triangular shaded regions F and G indicate the gravitational wave strain 
magnitude of extreme mass ratio inspirals (EMRIs) from binaries involving a supermassive black hole (SMBH) with (AI > 
10''A/q) and a smaller, compact partner. The upper region G bordered in black corresponds to signals produced by a IOMq 
black hole in orbit about a supermassive partner, while the grey-bordered region F corresponds to the signal from a O.SMq 
white dwarf in such a system. Vertical lines on the top horizontal border of the EMRI regions indicate the point at which the 
binary orbit shrinks to twice the gravitational radius of a SMBH with the indicated mass. The dark parallelogram E inset 
into the upper triangle corresponds to the gravitational strain magnitude of EMRIs involving the SMBH believed to lie at the 
center of Andromeda, for smaller partners with mass between 0.5Mq and IOMq. The dashed blue lines indicate the point at 
which inspiraling bodies double their frequency in less than one year (lower dashed line), one month (middle) or one day (upper 
dashed line). 



assume a background field of 1 mG is applied to fix the 
quantization axis, local fluctuations of 1 /LtG=10~^^ T per 
root Hz will cause phase errors of less than lO^^rad 
per root Hz, as required. This appears to be challeng- 
ing. We are primarily concerned with fluctuations on the 
timescale of T, and much less so with a constant back- 
ground. Magnetic field fluctuations on a time-scale of 
a second can be reduced to below 10~^'^T [60]. Thus 
in principle, magnetic fields can be kept sufficiently low, 
even for the optimized scenario. It remains challenging 
to achieve this over a tube length of kilometers rather 
than centimeters. 



2. Gravity gradient noise 

Gravity gradient noise [61] is caused by moving objects, 
such as vegetation, the sea, the atmosphere, transporta- 
tion, and most importantly, seismic activity. It cannot 
be shielded against, but it can be suppressed by choos- 
ing a location far away from from major perturbations. 
Unfortunately, relatively little is known about the exact 
magnitude of gravity gradient noise. To show how seri- 
ous this problem is, we shall present a simple estimate 
of seismic noise relevant to the basic interferometer. We 
will assume that it is installed with its highest point 2 
km below the surface. 

A convenient starting point for this discussion is the 
"new low noise model" for seismic noise on vertical com- 
ponent sensors presented in Fig. 1 of Ref. [62]. The 
dominant sources of noise are shown as being local at- 
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FIG. 9: Signals from extragalactic sources seen with an optimized AGIS. All sources are as indicated in Fig. 8. An optimized 
AGIS would be sufficiently sensitive to observe inspirals at the core of Andromeda. The dashed blue lines indicate the point at 
which inspiraling bodies double their frequency in less than one year (lower dashed line), one month (middle) or one day (upper 
dashed line). The orange region H denotes the signal produced by inspirals of WMq BH with a SOOM© to 10"* partner 
within ESO 243-49. Note that at distances approaching 1 Gpc, the effects of the cosmological redshift become significant, and 
this plot of source magnitude vs. frequency becomes inaccurate. Regions A-G are as previously defined. 
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mosphere noise up to about 2 mHz, hum due to high 
degree Earth normal modes near 3 mHz, Rayleigh waves 
near 10 mHz, and marine microseismics between 0.03 and 
2 Hz. According to Ref. [63], Sect. 6, "at quiet stations, 
much, and possibly most, of the vertical long-period seis- 
mic signal recorded in the period band 200 - 400 s has 
the dispersive properties of fundamental-branch Rayleigh 
waves." For estimating the transfer function from vertical 
displacement noise in this frequency band, we will thus 
model all seismic activity as Rayleigh waves. 

The vertical displacement due to a fundamental 
Rayleigh wave is given by [64] 



Tl 



Time 



FIG. 10: Space-time diagram of the BBB-interferometers. 1: 
Bloch oscillations in dual accelerated optical lattice; 2: Bragg beam 
splitter; 3: Bloch oscillations in quadruple lattice; 4; Dual Bragg 
beam splitter. 



= C (-O.Sbe-''''^ + 1.5e-'*'=(^^) cos k{ct 



(21) 



Here, z is the depth below the surface, x is the horizon- 
tal direction of propagation, k ^ 2tt/X the wavenumber, 
and c the speed of the wave. For an order of magnitude 
estimate, we use c — 205 m/s, q — 0.88 and s — 0.36 
as measured for the fundamental Rayleigh mode at the 
LIGO sites [611. The second term is dominant because 
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it is both larger and falls off more slowly with depth. 
Thus, we will assume the vertical displacement of matter 
is described by 



u = be ^'^^ cos k{ct - x) 



(22) 



where b is the surface vertical displacement. 

From this, wc can calculate the variation of the dis- 
tance to the earth's surface: The acceleration of free fall 
at the surface of an infinite sheet of thickness Z due to the 
sheet's mass is a = Tr^pGZ[\ — {Z/R)/'!t + . . .], where R is 
the radius of the sheet. The effect of the seismic wave is to 
replace Z z+u, so a = iT'^pG{Z+b)[l-{Z/R)/TT+. . .]. 
If the wavelength is infinite, this is going to have ex- 
act same effect on two atomic clouds and thus no effect 
on AGIS. To account for the finite wavelength, we use 
A ~ i?, and obtain a difference in accelerations of the top 
and bottom atomic clouds: 



at - aft = TTpGb- 



(23) 



where the local density p is about 1800 kg/m^ [61]. Fi- 
nally, Sa has to be divided by and by the baseline 
length L in order to obtain the resulting strain noise level 
in searching for GW signals, 



h = 



pGb 

2wc 



(24) 



Using the new low noise model [62] we obtain, for ex- 
ample, h = 2 X 10~^^/\/Hz at w = 27r X lOmHz, and 
h = 1.5 X 10"i^V\/Eb at w = 27r x lOOniHz. Both figures 
are between three and four orders of magnitude larger 
than the shot noise limit at 100 mHz (920 times as large) , 
and 10 mHz (14,000 times as large), respectively, than the 
shot noise (Tab. I and Fig. 3). The problem is even more 
severe for the optimized detector. 

This shows that gravity gradient noise is certainly one 
of the most important issues that need to be addressed 
to make AGIS viable. To suppress it, several strategies 
can be studies: 

• Use of more than two atomic clouds, whose signals 
are correlated in order to suppress seismic noise 
while retaining sensitivity to gravitational waves. 

• Additional suppression of gravity gradient noise 
can be obtained by seismic monitoring of the AGIS 
environment, but further study of the effects of 

seismically-induced gravitational gradients on an 
AGIS detector is needed. This can be expected 
to work best where it is needed most, for low fre- 
quency seismic waves. These waves have long cor- 
relation lengths and can be better modeled with a 
limited number of detectors. 

• Finally, setting up several AGIS can somewhat in- 
crease sensitivity through looking for correlations. 

However, the above very crude and simplified estimates 
show that the problem needs the attention of an expert. 



3. Influence of laser noise 

While a detailed analysis has yet to be performed, 
there are two requirements on the laser noise. The first 
is on low frequency fluctuations on the time scale of the 
pulse separation time T, the second is for fluctuations 
on the time scale of the duration t of the beam splitter 
pulse. 

High-frequency laser noise must meet the following re- 
quirement: An interferometer with a pulse separation 
time T at the shot noise limit requires the phase un- 
certainly per beam splitter be less than ^ \/ {n\/2rjT). 
If the beam splitting pulse takes a time a (half width), 
it samples the laser's phase noise spectral density over a 
bandwidth of roughly 1 / (27r(j) . The resulting phase noise 
spectral density of the beam would be 



\/2^ 



(25) 



if the two interferometers were independent. However, 
since a ^ L/c, this noise will mostly be common-mode 
to both interferometers. An estimate for the resulting 
requirement on the short term stability of the laser is 
thus 



v27ra tJ 



IQ-Vvlfe, 



(26) 



where the basic parameters (Tab. I) as well as cr = 100 ms 
were assumed. This is a mere -lOOdBc/v^Iz phase noise 
which is easily achievable at the state of the art of lasers 
[65]. Even the requirements for the optimized scenario 
(1.6 x IO-'^/a/Hz, or -135dBc/\/Hz) can be met. 

The low-frequency requirement arises since AGIS mea- 
sures the distance between two atoms by comparing it to 
the phase of a standing wave of a laser. In order to reach 
the atom shot noise limit, fluctuations bk in the effective 
wavevcctor must not fluctuate by more than the atom 
shot noise, n(Sk)L < l/^/2Tr}. The factor of n exists 
because the effective wavenumber is, to leading order, 
kes = nk. For the basic parameters in Tab. I, this leads 
to {Sk)/k < l/{nkLy/nj) « 4 x lO^^o^ rpj^jg corresponds 
to a ~ 12 pRz stability and must be maintained over the 
timescale of 2T. No such lasers exist at present: the best 
cavity stabilizations [66, 67] reach Hz-level stability. See 
Ref. [68] for prospects for a mHz laser. 

To illustrate this required level of stability, it is equiv- 
alent to the Doppler effect due to a c/ {nkL^2Trj) ^ 
1 X 10~^^ m/s velocity of the source. Thus, the position 
noise spectral density of the laser must not be more than 



2Tc 
nkL^ 



1.8 X 10-i°m/VIfe 



(27) 



on the time scale of 2 T. For comparison, the vibrational 

noise spectral density in the DUSEL underground facility 
was measured by Vuk Mandic (University of Minnesota) 
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to be around 10~^ m/v^Hz around 0.1 Hz, some three or- 
ders of magnitude short. 

To compare that to the corresponding requirement in 
a hght interferometer detector of the same sensitivity, we 
express x as a function of the sensitivity h in the low 
frequency hmit: x = chi.-pT^ijO^ ■ By comparison, a light 
interferometer detector requires a mirror position noise 
spectral density of roughly 5ligo = /i-^LiGO- Thus, the 
requirement for AGIS is less stringent than the one of 
light interferometers by a factor of 

5^ = ^~8xl0^f-^)' (28) 

a;iight Might V 2-77 Hz/ 

for Liight = 4 km and the basic scenario (the ratio is 
4 X 10^[a;/(27rHz)]^ for the optimized scenario). Thus, 
while position noise requirements for AGIS are extremely 
stringent, they are much less stringent than those for light 
interferometers. 

This points out a way in which the low frequency laser 
noise requirement can be alleviated: If it is possible to 
drive two AGIS sensors by the same laser, laser noise 
can be canceled. However, the beam splitting optics 
used must meet the vibration requirements just outlined. 
These requirements are much less stringent than corre- 
sponding requirements in LISA. Presenting a detailed 
scheme, however, is beyond the scope of the present work. 

F. Demonstrator Setup 

In order to demonstrate the critical technologies for 
GW detection, we are assembling a tabletop demonstra- 
tor in our lab. Our cesium setup is an atomic foun- 
tain, 1.5m tall [14, 16, 17], with a free evolution time 
T of 0.5 s. Inside our atomic fountain tube, we installed 
three layers of mu metal cylinders to suppress stray mag- 
netic fields. The atoms are trapped in a two-dimensional 
magneto-optical trap (2D MOT), generated with 2 W of 
laser power. This system should be capable of producing 
a flux of 10"'^^ atoms per second [49]. Atoms from the 
2D MOT are then transferred to a 3D MOT, and subse- 
quently launched by a moving optical molasses. The 3D 
MOT is about 2 cm in diameter and is estimated to con- 
tain on the order of 10^° atoms. The 2D MOT and 3D 
MOT chambers are separated by a differential pumping 
tube with a pressure ratio of 10"^, in order to maintain 
the high loading rate of atomic flux without increasing 
the pressure in the main 3D MOT chamber, where we 
are able to maintain a pressure of about 10"^" Torr. The 
lifetime of the atomic sample is thus on the order of a 
few seconds, which is crucial for interferometry. A tem- 
perature of less than 2 /xK has been achieved on a daily 
basis by using polarization gradient cooling and adia- 
batic cooling methods. Raman sideband cooling in an 
optical lattice [18] will be used to further cool the atoms 
to ~ 350 nK in the F = 3,mF = —3 quantum state. At 
this temperature, the atomic sample will expand to about 
1 cm over the course of the experiment. The atoms will 



be transferred to the mp = state by applying a small 
magnetic bias field and a lOW microwave sweep. A ve- 
locity selective Raman transition will reduce the vertical 
velocity width to 0.3 recoil velocities. The BBB beam 
splitter will be optimized for detection of gravitational 
waves. The aim is to achieve the high momentum split- 
ting required for AGIS. 

Our laser system for driving Bloch oscillations and 
Bragg diffraction will be based on a 6W Ti:sapphire 
laser [16, 57] (that we constructed from a Coherent 899 
laser without intracavity etalons) and a highly efficient 
arrangement of acousto-optical modulators. 

V. SUMMARY AND OUTLOOK 

We have studied candidate sources of gravitational 
waves for detection by an atomic gravitational wave inter- 
ferometric sensor, AGIS. We describe an optimized set of 
experimental parameters for reaching high sensitivity in 
a free- falling implementation of AGIS, where the atoms 
do not interact with the laser beams except during beam 
splitter operations. 

We consider binary inspirals to be the best candidates 
for detection by AGIS. A ground-based AGIS may be 
capable of detecting type la supernova precursors at up 
to 500 pc. An optimized detector should observe gravi- 
tational waves from at least one presently known binary 
system. AGIS may detect all neutron star-neutron star 
binary coalescences within 16 kpc, some up to a year in 
advance. It would also be sensitive to neutron star-black 
hole mergers within 50 kpc. AGIS may be able to observe 
the final decades of inspiraling black holes at the galac- 
tic core, and potentially the final years of an in-falling 
neutron star. AGIS would also be able to observe the fi- 
nal days of some extreme mass ratio inspirals within the 
Local Group, although it is uncertain that there are any 
supermassive black holes in the appropriate mass range. 
However, current estimates for the rate of occurrence of 
several classes of detectable events are spread over a dis- 
quietingly large range. 

Based on the shot noise limit, an optimized AGIS may 
be able to observe the late stage inspirals at the super- 
massive black hole in Andromeda, and to more general 
extreme mass ratio inspirals within the Local Superclus- 
ter, although we caution the reader that many other 
noise sources, such as Newtonian noise, may preclude 
this and remain to be studied. AGIS would then be use- 
ful for probing gravitational wave emissions from sources 
at non- negligible redshifts. 

However, several technical challenges need to be met 
in order to reach this: 

1. ultrahigh coherent momentum transfer to the 
atoms, 

2. high atom flux of lO^^/s, 

3. overcoming wavefront distortions in the laser 
beams. 



15 



4. high-power lasers will need to be developed. 

5. gravity gradient noise 

6. laser phase noise. 

Schemes for overcoming these disturbances need to be 
developed. These include a multi-arm version of AGIS, 
in which the effects of laser noise cancel, the use of mHz- 
linewidth lasers. Seismic noise seems to be about 3 or- 
ders of magnitude larger than the shot noise of even ba- 
sic AGIS, though strategies for reducing this might exist. 
For example, seismic monitoring and imdcrground oper- 
ation is currently being studied as a means of overcoming 
gravity gradient noise. This issue will definitely require 
the attention of experts. In addition, it has been pointed 
out [70] that wavcfront distortions in the laser beams to- 
gether with shot-to shot fluctuations in the location of 
the atomic cloud give rise to additional noise. 



Finally, we note that even our optimized scenario could 
be surpassed by more advanced atom-optics methods. 
For example, by holding the atoms in an optical lattice 
throughout the interferometer cycle, the pulse separation 
time can be extended beyond lis. 
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